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We review the criteria for separability and quantum entanglement, both in a bi- 
partite as well as a multipartite setting. We discuss Bell inequalities, entanglement 
witnesses, entropic inequalities, bound entanglement and several features of mul- 
tipartite entanglement. We indicate how these criteria bear on the experimental 
detection of quantum entanglement. 
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The phenomenon of quantum entanglement lies at the heart of quantum me- 
chanics. And what lies at the heart of quantum mechanics, may lie at the heart 
of a future technology. It is not surprising then that over the last 5 years a 
theory of quantum entanglement has started to emerge that tries to capture, 
quantify and assess the power of quantum entanglement. 

First, it was by the protocol of quantum teleportation [BBC + 93] that quantum 
entanglement was introduced as a resource in quantum communication: it has 
become a rule of quantum communication law that one bit of entanglement 
(1 ebit) enables 1 unknown qubit to be sent by means of 2 classical bits. 

But it has been realized over the last year that quantum entanglement is 
not only a fundamental resource in quantum communication, but can also 
be viewed as a resource in quantum computation. Gottesman and Chuang 
[GC99] have shown that it is possible to perform universal quantum compu- 
tation, by starting with three-party entangled GHZ states and subsequently 
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performing single qubit operations and measurements in the Bell basis. In the 
linear optics quantum computation proposal by Knill, Laflamme and Milburn 
[KLM] the quantum gate that lies beyond the capabilities of linear optics, 
can in fact be implemented by the creation of a multipartite entangled state. 
Quantum entanglement also lies at the core of the quantum computation pro- 
posal by Rausschendorf and Briegel [RB]. In this proposal the authors show 
that universal quantum computation is possible by means of a series of single 
qubit measurements that are performed on an initial state which is a certain 
highly entangled 'cluster' state. The entangled state functions as a substrate 
on which the quantum computation takes place. 

In this article we review the progress that has been made in establishing one 
of the cornerstones of the theory of quantum entanglement, namely the de- 
velopment of criteria for entanglement and separability, both in the bipartite 
as well as the multipartite setting. In the last section of this paper we will 
consider how entanglement witnesses can be used in deciding by experiment 
whether a quantum state is entangled. We will not discuss the topic of entan- 
glement measures, which can be viewed as a subject complementary to the 
one which we consider in this review article. We would like to refer the reader 
to Ref. [HHHOO] for a more comprehensive overview on bipartite quantum 
entanglement. 

We will write X, Y and Z for the three Pauli matrices. A positive semi definite 
operator A with nonnegative eigenvalues is denoted as A > 0. A n-dimensional 
Hilbert space is denoted as H n , and operators on this space (n x n matrices) 
G B(H n ). Furthermore the class of quantum operations which are constructed 
by Local Operations supplemented by Classical Communication is sometimes 
abbreviated as LOCC. 



1.1 What is Quantum Entanglement 

A bipartite pure quantum state \ip) G Ha <8> He is called entangled when it 
cannot be written as = \ip a) <8> \^b) for some \iPa) G Ha and G He- 
A mixed state or density matrix p, which is a positive semidefinite operator 
on the space Ha ®Hb, is called entangled when it cannot be written in the 
following form 

p = Eftl^X^I®l^f>^fl. (!) 



for some set of states \ipf) G Ha, £ He and pi > 0. If the density matrix 
p can be written in the form of Eq. (1), then the density matrix p is called 
separable. 
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The distinction between separable states and entangled states has an oper- 
ational meaning in the following sense. A source (or black box) produces a 
mixed state p; the mixedness could come when aside from systems A and B 
there are additional degrees of freedom to which we have no access. If p is 
an entangled density matrix, then some coherent interaction must have taken 
place between A and B. If p is separable, then no guarantee exists of whether 
the interaction in the black box was coherent or not. 

Consider a multipartite system with parties labeled by A±, . . . A n . The density 
matrix is called separable when no entanglement exists between the parties, 
i.e. p = J2i Pilipf 1 } (ipf 1 ]® ■ ■ •®Vl } t n ){' t l } t n \- F° r the various degrees and forms of 
quantum entanglement that can exist among parties, we will consider specific 
classes of states in Section 3. 

A note of caution about how to interpret the state of a physical system in 
terms of quantum entanglement may be in place here. The previous standard 
definitions of quantum entanglement tacitly assume that (1) every state in 
the bi- or multipartite Hilbert space is in principle available as a physical 
state and (2) local (involving single tensor factors) as well as global quantum 
operations, measurements and unitary transformations, can be performed on 
the Hilbert space. In this respect the wavefunction of two identical bosons 
^/(xi,X2) = ip{xi) < 8>ip( x 2) + '4 , { x 2) <S>ip(xi) cannot be called entangled, since it 
falls short of these criteria. Understandably, when considering more complex 
physical systems, the dividing line between what is entangled and what is 
not entangled, may become somewhat fuzzy. The guideline in deciding these 
matters, I believe, should be the question: "Do we have an operational form 
of quantum entanglement? What resource does the particular state constitute 
in quantum communication and computation?" In Ref. [SCK + ] for example 
the authors consider the entanglement that can exist in 2-fermion systems. 



2 Bipartite Criteria 

2.1 Bell Inequalities 

Historically one can say that the first separability criterion was formulated 
by John Bell [Bel64]. Bell's intention however was not to establish a separa- 
bility criterion, but to evaluate the power of local hidden variable theories in 
describing local measurement outcomes on quantum mechanical states. 

His inequality, and similar inequalities such as the CHSH inequality [CHSH69] 
found later, is obeyed by any local hidden variable theory, whereas the cor- 
relations in measurement outcomes on, for example, the singlet state = 
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^(|01) — 1 10)) violate the inequality. Furthermore the outcomes of local mea- 
surements on any separable density matrix can be simulated by a local hidden 
variable theory. This can be easily understood from Figure 1 which gives an 
idea of the workings of a local hidden variable theory. For every pure entangled 
state there exists a Bell inequality that is violated [Per93] and therefore there 
exists a series of measurements and outcomes through which we can ascertain 
that our state is entangled. 



Fig. 1. Local Hidden Variable Theories: Alice and Bob each have an arbitrarily 
powerful machine in their lab which takes as input the description of their local 
measurements Ma and M b and a description of the state p on which the measure- 
ment will take place. Inside their machine may be a random shared bit string r of 
arbitrary length. The output of the machines is supposed to statistically simulate 
the outputs of the real measurements Ma and Mb that were performed on the 
state p, in the sense that the joint and marginal probabilities for various outcomes 
and choices of measurements are identical to the real measurement on p. 

The weakness of Bell inequalities as criteria for entanglement or separability 
lies in the fact that it is not known whether violations exist for many entangled 
mixed states. For example, it has been shown that for a special class of mixed 
states, so called PPT bound entangled states (see Section 2.2), all CHSH- 
inequalities are obeyed [WW00]. If we loosen the rules of the game and allow 
preprocessing of our state p or many copies of our state p® n by means of 
LOCC, then a much larger class of states p will violate a Bell inequality. 
We demand that only local operations and classical communication enter in 
this game, since these are the operations which cannot increase the quantum 
entanglement in a state. This class of operations, crucial in the theory of 
quantum entanglement, is graphically depicted in Fig. 2. All density matrices 
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Fig. 2. Local Operations and Classical Communication: at each round Alice's (Bob's) 
local actions may depend on Bob's (Alice's) previous actions and outcomes. 

which are distillable (see Section 2.3) will then violate a Bell inequality in this 
manner. 



2.2 Entanglement Witnesses and Positive Linear Maps 

The framework of Bell inequalities fits in a larger scheme of entanglement 
witnesses. In fact, each Bell inequality can be viewed as a particular example 
of an entanglement witness [TerOO]. A prime example is the operator form 
of the CHSH inequality: the Bell-CHSH operator (a, a', b, b' are unit-vectors) 
reads 

B = a ■ a <g> (b + V) ■ a + a> ■ a <g> (b - b') ■ a. (2) 

The expectation value of B with respect to all separable states 

Tr£p sep <2, (3) 

whereas Tr Bp can exceed this value for an entangled state p. The operator 
21 — B is an example of an entanglement witness. Even though there does 
not necessarily exist a Bell inequality for every entangled state (for certain 
Werner states, for example, there is no single copy violation of a Bell inequality 
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[Wer89]), there does exist a witness for every entangled state. This is the 
content of the following theorem: 

Theorem 1 (Horodecki [HHH96]) A density matrix p on Ha ®7~Lb is entan- 
gled if and only if there exists a Hermitian matrix H = H\ an entanglement 
witness, such that 

TrHp<0, (4) 

and for all separable states p sep , 

TrHp sep >0. (5) 

Thus by a measurement of the entanglement witness observable H we will be 
able to decide whether a particular state is entangled, since when we find a 
negative expectation value for H, we must conclude that the state cannot be 
separable. 

There exists a direct relation between entanglement witnesses and positive lin- 
ear maps which are not completely positive [Jam72]. A linear map C: B(H n ) — > 
B(H m ) is called positive, when it maps all X > onto C(X) > 0. The map £ 
is completely positive if and only if l n ®£ is a positive map. The most famous 
and physically relevant example of a positive map is matrix transposition T 
in a chosen basis. The map T is not completely positive, as can be illustrated 
by applying it on half of a (unnormalized) maximally entangled state: 

(i®r)[|oo> + |n»«oo| + (ii|] = 

|00)(00| + |11)(11| + |01)(10| + |10)(01|. (6) 

The resulting operator has an eigenvector |01) — 1 10) with negative eigenvalue 
— 1. It was noted by Peres [Per96] that applying 1 <8>T on a separable density 
matrix always gives another density matrix 

i 

E^I^)(^l®I^D(^ri>o, (7) 

i 

and therefore the condition (1 <g> T)(p) > 0, sometimes called the Peres- 
Horodecki criterion, constitutes a separability criterion. 

The relation between entanglement witnesses and positive linear map is the 
following. We take a maximally entangled state in H n <8> H n , for example the 
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state |^+) = £™=i The Hermitian operator H e B(Ti n ®Ti m ) defined by 

H={l®C){\^ + )^ + \), (8) 

has the property of Eq. (5) if and only if C: B{7i n ) — > B(7i m ) is a positive 
map. Furthermore, H is an entanglement witness, as in Eq. (4), if and only if 
C is not a completely positive map. 

The theory of positive maps which are not completely positive has not been 
completely developed. What is known is that in spaces such as H.2 <8> 7i,2 an d 
Ti.2 <8> 7~(-3, all positive maps C relate to the matrix transposition map T, i.e. 
they can all be written as 

C = S 1 + S 2 o T, (9) 

where <Si and S 2 are completely positive maps and T is matrix transposition 
in any chosen basis [Wor76]. This implies that for these small dimensions, 
Ti-2 <E> Ti-2 or 7i 2 <E> H3, the entanglement witnesses are of the form 

H = P + (1®T)(Q), (10) 

where Q > 0, P > and T is matrix transposition in any chosen basis. 

In higher dimensions the situation is more involved. There do exist positive 
maps which are not decomposable, meaning that they are not of the form of Eq. 
(9). A consequence is that in higher dimensions, there exist entangled states 
which satisfy the Peres- Horodecki criterion, i.e. (1 <8> T)(p) > where T is 
matrix transposition in any basis. These are the bound entangled density ma- 
trices with PPT (Positive Partial Transposition) (see the next Section). The 
first example of an indecomposable positive map in TL3 <E> TC3 was found by 
M.-D. Choi [Cho75]. The first method of constructing some indecomposable 
entanglement witnesses in arbitrary high dimensions was presented in Ref. 
[Ter]. In Refs. [LKHC] and [LKCH] this construction was generalized with the 
following consequences. It was shown in Ref. [LKHC] that every indecompos- 
able entanglement witness is of the form 

H = P + (l®T)(Q)-el, (11) 

where, in order to ensure that H has the property of Eq. (5), we have 

0<e< inf {^ A ^ B \(P+(1®T)(Q))\^ A ^ B ). (12) 

WA,WB 
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Here P > 0, Q > and they are such that Tr P5 = and TrQ (1 <g> T)(S) = 
for an 'edge' state 5. An edge state 5 is a bound entangled PPT state 
which has the property that for all e > and all product states I^a^b), 
5 — e\ipA, i J B){i J A, i/>b\ is not positive or does not have PPT. As such, the edge 
states are on the boundary of the closed set of entangled PTT states. The 
entanglement witness H in Eq. (11) detects the entanglement of the edge 
state 5. 

The entanglement of bound entangled PPT states p in the interior of the 
set of PPT states is also detected by these witnesses [LKHC]. By choosing 
Tr(P + (l<S>T)(Q))p = we ensure that the indecomposable witness H in Eq. 
(11) detects the entanglement in p, i.e. TiHp < — e < 0. 

Thus, given an edge state 5, its entanglement witness can be determined. A 
complete characterization of these edge states is still an open question; in Ref. 
[LKHC] they are shown to be based on pairs of subspaces Tii and H.2 such 
that (1) for every product state IV'AjV'b) £ 'Hi, I^a^b) & ^2 and (2) the 
rank of Tr i P-^ 1 is equal to the rank of TrjP-^ 2 for % — A, B, where P is the 
projector onto the subspace. For the choice Hi = H2, an example of such a 
subspace (containing, in this case, no product vectors) is the space orthogonal 
to an unextendible product basis [BDM + 99]. It would be very interesting to 
find a method for constructing such subspaces in general. 

Indecomposable positive maps are highly nontrivial objects. As was noted by 
M.-D. Choi they provide special counterexamples to Hilbert's 17th problem. 
In Appendix A we will review this connection. 



2.3 Operational criteria: LOCC and distillability 

The resource view of quantum entanglement emerged with the discovery of 
quantum teleportation [BBC + 93]. A natural next question, which was asked 
and partially answered by Bennett et al. in Ref. [BDSW96] is, in what sense 
does mixed state entanglement enables quantum data transmission? This line 
of study, which is still ongoing, has led to operational criteria for 'useful' 
entanglement. We can classify entangled states in terms of their distillability. 
Distillation of a mixed state p is a process which is implemented by LOCC 
on a large set of copies of the state p® n . The process outputs a smaller set of 
states a® k on a space of dimension 2 k x 2 k such that (^" 0fc | a® k \^~® k ) -> 1 
when n — > 00. Here \^~) is the singlet state. The asymptotic fraction - is 
called the distillable entanglement D of the density matrix p. It has been 
shown that all entangled density matrices in 7i 2 ®7i2 are distillable [HHH97]. 
In fact, a generalization of this result exists: all states which violate the so 
called reduction criterion are known to be distillable [HH99]. The reduction 
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criterion is violated for a density matrix p when either 



1®Pb~P^0 or p A ®l- P^O. (13) 

It is noteworthy that satisfaction of the reduction criterion is identical [HH99] 
to positivity under partial application of the decomposable positive map 
C(X) = lTrX — X, since (1 ® C)(p) = Pa <8> 1 — p- Since the map C is 
decomposable, Eq. (9), it follows that every state which satisfies the Peres- 
Horodecki criterion will also satisfy the reduction criterion. 

The class of density matrices which are not distillable are called bound en- 
tangled density matrices. At least one group of density matrices exists for 
which it is possible to rigorously prove that they are bound entangled den- 
sity matrices. These are entangled density matrices p which do not violate 
the Peres-Horodecki criterion. This criterion is important in the theory of 
quantum entanglement, since it has been shown [HHH98] that nonviolation 
of the Peres-Horodecki criterion is preserved under LOCC. This is not hard 
to see, when we consider a larger class of quantum operations, the separable 
superoperators, of which the LOCC actions are a subset. A separable super- 
operator or measurement [Rai99] has operation elements that are of separable 
form {Ai <g> B u Ei A A i ® B \ B i = !}• Given is a density matrix p such that 
(1 <S> T)(p) > 0. For every i, we can write 

(1 <g> T)[A <g> B iP A\ ® B\\ = (Ai <g> B*)(l ® T)(p)(4 <g> Bj) > 0, (14) 

under the premise that p satisfies the Peres-Horodecki criterion. Thus under 
the action of a separable superoperator or measurement, PPTness is preserved. 
Since the output of a distillation process (arbitrarily good approximations to 
singlet states) does violate the Peres-Horodecki criterion, it must follow that 
entangled states with the PPT property cannot be distilled. The first bound 
entangled states were found by P. Horodecki [Hor97]. In Refs. [BDM + 99] and 
[DMS + ] constructions were given for classes of bound entangled states based 
on unextendible product bases. 

A second class of states has been conjectured to be nondistillable 
[DSS + 00,DCLB00]. Examples of these states are particular Werner states in 
7~t n ® H n , of the form 

Pa = A(n2 _\ ) _ 1 (Al - (A + 1)(1 ® T)(|*+)(*+|)) , (15) 

where \^ + ) = ^= Y%=i \hi)- F° r an finite A > 0, these Werner states violate 
the Peres-Horodecki criterion. The states are conjectured to be nondistillable 
for A G [^2,00). 
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The structure among the set of (conjectured) nondistillable states is by it- 
self fairly complex. It has been shown recently that the tensorproduct of two 
nondistillable states, one of the PPT kind and one nondistillable Werner state, 
Eq. (15), in H 3 ® H3, can be a distillable state [SSTa]. Furthermore, it is not 
clear what classes of bound entangled states are asymptotically interconvert- 
ible by local actions and classical communication. In order to investigate this 
'fine-structure' we may need to look for positive maps C for which "positivity 
under partial application of £" is preserved under LOCC. 

2.4 Functional Separability Criteria 

Even though entanglement witnesses completely characterize the set of sepa- 
rable states, they do not provide a simple computational method (except for 
the Peres-Horodecki condition) for deciding whether a density matrix is en- 
tangled. It is desirable to have alternative or additional criteria that may help 
us in deciding this and characterizing separable versus entangled states. The 
criteria that we will discuss in this Section are all obeyed by separable density 
matrices, so that in case of a violation we know that p is entangled. 

The first criterion is a combination of the Peres-Horodecki criterion and a 
check on the rank of the density matrix. When we find that p has PPT, p may 
either be separable or bound entangled. In Ref. [HLVCOO] it was proved that a 
density matrix p G B{7i n ®l-L in ) with PPT and a rank which is smaller than or 
equal to max(n, to) is separable. For PPT density matrices for which the sum 
of the rank of p and (1 ® T)(p) is smaller than or equal to 2mn — m — n + 2, 
the authors in Ref. [HLVCOO] provide an algorithm for checking whether p is 
separable. 

It can be shown that from a non-violation of the reduction criterion several 
other criteria can be derived: 

Lemma 2 For all separable states the reduction criterion is not violated and 
this implies that if p is separable, 

S a {pA) < S a (p) and S a {p B ) < S a {p). (16) 

for a — 0, a G [1, 2] and a = 00 where S a (p) for < a < 00 (a ^ 1) is the 
quantum Renyi entropy 

S a (p) = -^—\ogTrp a . (17) 
I — a 

For a = 0, we have S (p) = log R(p) where R(p) is the rank of p, lim^x S a = 
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S(p) where S is the von Neumann entropy S(p) = — Trplogp and for a = oo 
we have = — log | \p\ \ . 

The case a = oo was proved in Ref. [HH99] and the case a = was proved in 
Ref. [HSTT]. The cases a — > 1 and a = 2 can be derived from the reduction 
criterion [HH]. For a — > 1, we use the operator- mo notoniticity [And94] of the 
log function to infer that 

log p A ® 1 > log p. (18) 

Now we use that when X > 0, Tr pX > for all p > 0. Thus, multiplying with 
p on both sides and subsequently taking the trace gives the desired entropy 
inequality for a — * 1. If Pa®1 > P, then (p^£g>l) 5 > p" 5 for 5 e (0, 1], since t 5 is 
operator- monotone for 5 in this interval [And94] . If we multiply the inequality 
with 5 by p on both sides and trace, we get 

Tr A p A +s > Trp 1+6 , (19) 

which, after taking logarithms on both sides, proves the result for a G (1, 2]. 
In Ref. [CAG99] the conditional quantum operator was defined 

p A \ B = exp [log p- log 1 A Ops]- (20) 

With this definition the conditional entropy S(A\B) = — Trplogp^is is the 
difference between the total entropy of the state S(p) and the local entropy 
S(pb)- Lemma 2 states that for all states obeying the reduction criterion (sep- 
arable states and at least all nondistillable states) such a conditional entropy 
(and similarly some a-entropic extensions) is nonnegative. 

Nielsen and Kempe [NK] recently found a different separability criterion 

Lemma 3 [NK] For all separable density matrices p 

\ PA y \ p and \ PB y X p , (21) 

where X a is the ordered vector of eigenvalues of the density matrix a. Here the 
symbol >~ means majorization, i.e. A >- p, when Y$=i \ > Z)£=i A*i f or a ^ k- 

Similar to the entropic criteria given above, this majorization criterion af- 
firms the intuition that separable density matrices are globally at least as 
mixed as locally. A new corollary of the majorization criterion are the en- 
tropic inequalities in Eq. (16) for a G [0, 1]. This follows from the fact that 

— # 

the quantum Renyi entropy S a (p) is a concave function of the probabilities A p 
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for all a G [0, 1]. In Ref. [Well 78] it is claimed, in fact, that the majorization 
condition (Uhlmann's relation) implies the entropic inequalities for all a > 0. 

The two criteria, the reduction criterion of Eq. (13) and Lemma 3 are strik- 
ingly similar. There exist states however for which the reduction criterion is 
violated whereas the majorization criterion is satisfied; this is the example of 
a 2 qubit entangled state in Ref. [NK]. The reduction criterion is violated for 
this state, since in 7i 2 ® 7i 2 the reduction criterion is equivalent to the Peres- 
Horodecki criterion, and all entangled 2-qubit states violate this criterion. The 
conjectured nondistillable Werner states, Eq. (15), obey both the reduction as 
well as the majorization criterion. It is possible, but unproven, that any state 
which satisfies the reduction criterion also satisfies the majorization criterion. 

Unfortunately, non-violation of the reduction criterion and non-violation of 
the majorization criterion are not properties that are preserved under local 
actions and classical communications. To take an example, consider the fol- 
lowing density matrix a = 1a/ d® \^! + ){^ + \ <8> 1b /d where 1a/ d (1b/ d) is a 
density matrix for Alice (Bob), is a maximally entangled state in TC n ®H n 
and d is large. The state P + = violates the reduction criterion 

1/n - P + t 0, (22) 

We can always choose d > n large enough such that both the reduction criteria 
are satisfied for a: 

1 A <S> [1a <S> ls/n - P+/d] <g> 1 B > 0, and 

1 A ® [l A /n®l B -P + /d\®l B >0. (23) 

By the local action of tracing over the register with l A /d and ls/d, we obtain 
the state P + which does violate the reduction criterion. Thus there exists 
states which initially do not violate the reduction criterion, but nonetheless 
are distillable. The same example can serve to show that the majorization 
criterion can be violated only after some local action. 

In Figure 2.4 we have sketched an overview of the known relations between 
the various bipartite separability criteria. 



3 Multipartite Entanglement 

Every time we consider a multipartite system in a bipartite fashion, that is, we 
split the set of parties in two subsets, and consider the entanglement between 
the subsets, we can apply the criteria that we have listed in the previous 
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Fig. 3. Relations between various bipartite separability criteria. It is not known 
whether the green (Reduction Criterion) and light blue (Peres-Horodecki) ellipses 
are contained in the dark blue (Majorization) ellipse. In Ti.20Ti.2 all colored ellipses, 
except the dark blue one, collapse onto each other and mark the separation of the 
set of separable states from the entangled states. 

Section. In this Section we will focus on features of multipartite quantum 
entanglement which are special to multipartite quantum entanglement. 

3.1 Violations of Local Realism and Bell inequalities 

The Greenberger-Horne-Zeilinger state \GHZ) = ^=(|000) — 1 1 1 1} ) is an ex- 
ample of a three party state which violates the predictions of local realism 
[Mer90b]. The GHZ-state is an eigenvector of operators X <g>Y <g>Y , Y (g> X <g>Y , 
Y ®Y ® X with eigenvalues +1 and X ® X ® X with eigenvalue —1. These 
first three operators form the generators of an abelian group G which contains 
X <S> X <g) X. For these four operators, from a measurement of X or Y on 
two of the spins we can deduce the outcome of the third since \GHZ) is an 
eigenstate. Therefore, according to local realism, we may assign values to the 
local operators Xj and Yi according to some function / 

f:X t ^ {-!,+!} f:Yi^ {-1,+1}, (24) 

where X { acts on the ith particle, obeying the eigenequation constraints. The 
function / gives rise to a function h: g E G — > { — 1, +1}, i.e. h(X <g> Y <8> Y) = 
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f (X) f (Y) f (Y) . A violation of local realism occurs when it is impossible to 
construct a local function / (consistent with the eigenvalue equations for the 
generators) such that h is a group homomorphism, i.e. h(g 1 og 2 ) = fo(<?i)o/i(<? 2 ) 
for all <?i,<72 £ G. In the case of the GHZ-state, the violation comes about 
by observing that the local assignments (the function f) for the generators 
X ®Y ®Y, Y <g> X <g> X and X ® X ®Y always give h(X <g> X <g> X) = 1 
whereas the GHZ-state has eigenvalue —1 with respect to X <g> X <g> X. 

In Ref. [DP97] violations of local realism were found for a special class of 
multipartite entangled states. These are states used in quantum error correc- 
tion codes based on the stabilizer formalism. The states are by definition the 
eigenvectors of an Abelian group made from tensor products of Pauli matrices 
and 1. 

These violations do not yet present us with a separability criterion. The only 
claim is that if the outcomes of local measurements of X and Y were to corre- 
spond exactly to, say, the outcomes on the GHZ-state, then we may conclude 
that these measurement outcomes cannot be described by a local hidden vari- 
able theory. However, to establish a full separability criterion, we would need 
to analyze what ranges of outcomes could still be reproduced by a separable 
state or local hidden variable theory and what outcomes cannot. 

The n-qubit Bell-Klyshko operator [GBP98] or Mermin's inequality [Mer90a] 
in operator form (see for example Ref. [WWOO]) do constitute a separability 
criterion in this way. The Bell-Klyshko operator can be defined recursively 

1 1 

B n = £ n _i <g> -(a n ■ a + ct n ■ a) + B' n _ x <g> -(a n ■ a - a! n ■ a), (25) 

where B' n _ l = B n -i(a\ <-> a[, ai <-> ct 2 , ■ ■ ■ , a n -i ^n-i) an d B2 is the Bell- 
CHSH operator in Eq. (2). Its expectation value for all separable states is 
bounded as Tr B n p sep < 2, whereas a value as large as 2^ n+l ^ 2 can be obtained 
for the cat state \GHZ n ) = ^(|0® n ) + |l 0ri )). As in the bipartite case, these 
operators are examples of entanglement witnesses. 

3.2 Entanglement Witnesses and Linear Maps 

In Ref. [HHH] the notion of entanglement witnesses was formally extended 
to the domain of multipartite quantum systems. The multipartite analog of 
Theorem 1 separates the multipartite separable states from any entangled 
state; for every state p, not of the form p = Y^iPi^f 1 ) (ipf 1 \ <8> • • • <8> \i ) t n )('4 , t n \i 
there exists a Hermitian operator H, the entanglement witness, such that 

TrH Psep > 0, (26) 
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for all separable states p and Tr Hp < 0. Interestingly, these multipartite 
witnesses can be shown [HHH] to relate to linear maps C: B(7i 2 ® H? ® . . . <g) 
TC n ) — > B{Ti}) with the property that for all product vectors |x 2 , . . . , 

£(|x 2 , . . . ,x n )(x 2 , • • ■ ,x n \) > 0. (27) 

The linear map C is thus not necessarily positive on B(H 2 ®'H Z ®. . .<8>7i n ); this 
is an important difference with the bipartite case. The 1-1 relation between a 
witness H with the property of Eq. (26) and the map C with the property of 
Eq. (27) is the following: 

H = (l 1 ®tf)(\* + ){* + \), (28) 

where is a maximally entangled state in 7i l <g> Here the Hermitian 
conjugate of C, O is defined by the relation Tr A*£(B) = Tr tf(A*)B. 

In Ref. [JB00] a family of multiqubit entanglement witnesses was presented. 
For an n-partite qubitsystem the authors introduce an averaging observable 
a — - Y%=i a i where is an observable acting on the ith factor in H} ® Ti. 2 ® 
. . . <E> H n . The operators bounded in their norm, ||aj|| < 1 (here ||.|| is 

the standard operator norm). The operator a could, for example, be a sum of 
Pauli operators Z it measuring a mean magnetic field in the z-direction. It is 
proved in Ref. [JB00] that for separable density matrices p the expectation of 
the commutator 

|Trp[a,c]|<4=, (29) 



for any averaging observable a and c — <y with ||c|| < 1. An expectation value 
such as Eq. (29) can appear in a perturbative expansion in linear response 
theory (see for example Section IV in Ref. [TD00]). The operator c will then 
be a nonlocal time-dependent observable, the operator a some local mean 
field observable and p can be the equilibrium state of a physical system. For 
entangled states the expectation value of the commutator in Eq. (29) can be of 
order 1 as was shown in Ref. [ JB00] . A consequence is that for such entangled 
states the region of validity of the perturbation theory is much smaller than for 
separable states. Let us translate the result in the language of entanglement 
witnesses. The observable 

H = -=l-i[a,c], (30) 

is a witness in the sense that Eq. (26) holds. This witness can detect the 
entanglement in superpositions of macroscopically distinct states, such as the 
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cat \GHZ n ) state. For example (see Ref. [JBOO]) when we choose the operators 
at = and c = i[|O® n )(l 0n | - |l® n )(O 0n |], we obtain a witness which has 

(GHZ n \ H\GHZ n ) = -%=-!, (31) 



which is negative when n > 2. 



3.3 Incomparable Forms of Pure State Entanglement 

The possibility for extraction of bipartite pure state entanglement from mixed 
state entanglement by LOCC is captured by the notion of distillation, which 
we discussed in Section 2.3. The definition of bipartite distillation does not de- 
pend on the form, -maximally or partially entangled states-, of the final pure 
state entanglement, since the Asymptotic Interconversion Theorem for bipar- 
tite entanglement [BBPS96] says that all bipartite pure state entanglement is 
interconvertible by LOCC in the asymptotic limit (when we have many copies 
of a state). For multipartite entanglement no such theorem exists. The explo- 
ration of interconvertibility of multipartite entanglement was initiated in Ref. 
[BPR+00]. It was found for example that 3 EPR pairs (|00) + \ll)) m shared 
among three parties were not exactly convertible by LOCC to 2 GHZ states 
(| 000) + | 111)) 02 . In Ref. [LPSW] this result was considerably strengthened by 
showing that 3n EPR pairs are not interconvertible to 2n GHZ states even in 
the asymptotic limit n — ■> oo. As such these states form the building blocks of 
an MREGS, a Minimal Reversible Entanglement Generating Set [BPR + 00], 
with which all tripartite entanglement can reversibly be created. It is an open 
question whether a third type of state, the H^-state |001) + |010) + 1 100) 
should be added to the tripartite MREGS, in other words whether the Ha- 
state is asymptotically interconvertible to a supply of EPR pairs and GHZ 
states (see Refs. [DVC00] and [GPV] for indications that this may not be the 
case). These results show that intrinsically multipartite forms of pure state 
entanglement exist. 



3.4 Bound Entanglement 

One of the best illustrations of the phenomenon of intrinsic mixed state 'mul- 
tipartiteness' was given in Ref. [BDM + 99]. It is an example of a tripartite 
mixed state in 7^2 ® ® Ti.2 which is separable over all bipartite cuts of the 
three parties, while at the same time the state is entangled. Let \vi) = |000), 
\v 2 ) = I - +1)> 3 > = 1 + 1") and \v 4 ) = |1 - +), where |±) = -L(|0) ± |1)). 
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The state is 

4 

Pshifts = 1 - ^2\vi)(vi\. (32) 
i=i 

No product state state exists in the range of pshifts, since the vectors 
{\vi), \v2), 1^3), 1^4)} form an unextendible product basis [BDM + 99]. At the 
same time, the product basis can be completed with vectors which are separa- 
ble over a bipartite cut, which results in p being separable over this cut. The 
density matrix pshifts m Eq. (32) is also an example of a bound entangled state 
in the multipartite setting; if entanglement could be distilled from /ofhifts by 
local actions and classical communications of the three parties, then entangle- 
ment would be created over some bipartite cut, which is forbidden since p is 
separable over all cuts. 

The phenomenon of multipartite bound entanglement is more general than 
this. It can be argued [Smo] that any multipartite density matrix p for parties 
Ai, . . . , Ak is bound entangled when for all pairs of parties (Ai, Aj) there exists 
a cut (a bipartition) where Ai and Aj are in different sets of parties, such that 
p is either separable or PPT over this cut. This follows from the fact that these 
separability and PTT properties are preserved under LOCC and the fact that 
any multipartite pure entangled state is entangled over some bipartition. In 
Ref. [DCOO] and Ref. [SSTb] multiqubit examples of such bound entangled 
states were presented. Moreover in Ref. [SSTb] it was shown that two 4-party 
bound entangled states both distributed among 5 parties, can be distillable. 
The state p ACBD is a mixture of Bell states; with probability 1/4, A and C 
share one of the 4 Bell states and B and D share the same Bell state. The state 
is symmetric under permutation of parties. Furthermore it is separable over 
any bipartition into (2,2), but it is entangled for any bipartition into (3, 1). 
The other state p ABCE is identical, except that now it is shared among A, 
B, C and E. The essential feature of these two states taken together, is that 
there exists no bipartition over which p ACBD ® p AB CE j g separable S uch that 
the parties D and E belong to different sets. This fact makes it possible for 
entanglement to be distilled between D and E and the authors of Ref. [SSTb] 
show that 1 singlet (1 ebit) can be obtained from the two bound entangled 
states. 



4 Experimental Issues: Detecting Quantum Entanglement 

In this Section we consider how the criteria that we have discussed in the 
previous Sections enable us to decide by physical experiment whether the 
physical state of a given quantum state is entangled. The capacity to build 
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certain entangled states is one of the basic requirements for making a quantum 
computer and is often used as a benchmark test for the amount of control and 
coherence in a particular quantum system, see for example the creation of the 
cat state |0® 7 ) + |1® 7 ) in the NMR experiment in Ref. [KLMTOO]. 

It is desirable that the verification of the entanglement take place with a min- 
imal number of measurements and operations. The first, but inefficient option 
would be to perform full quantum tomography on the state p, i.e. determine 
all the matrix elements py, after which we may analyze the state on paper 
using the various separability criteria. If no prior knowledge exists about the 
state, then it appears that there is no shortcut to such a quantum tomography 
experiment. In the more common situation in which we expect to have created 
a certain state p, more efficient methods exist. A traditional method for de- 
tecting entanglement (employed for example in Ref. [BPD + 99]) is to test for a 
violation of a Bell type inequality. We have indicated in this review that these 
tests are part of a larger framework of entanglement witnesses which exists 
both in the bipartite as well as in the multipartite setting. In Ref. [LKCH] the 
notion of an optimal entanglement witness was introduced; an entanglement 
witness H is optimal if there exists no other witness H' which detects the same 
entanglement (the same states) as H and more. These optimal witnesses are 
the ones that will be useful in detecting quantum entanglement. We will need 
the following additional definition: 

Definition 4 A p-optimal witness for an entangled state p is an optimal 
entanglement witness according to the definition in Ref. [LKCH] and among 
the optimal normalized witnesses, it is the best in detecting the entanglement 
of p, i.e. 

Tr H*p = min Tr Hp, (33) 
optimal H 



with H normalized as 

TiH = 1. (34) 



In the next Sections we will consider the optimal ^-witness for pure states 
\ip) and small systems, for mixed states in larger systems and for multipar- 
tite entangled states. In Section 4.4 we consider how well these entanglement 
witnesses detect entanglement in the vicinity of the desired entangled state. 
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4-1 (Multipartite) pure states, small bipartite mixed states 

Let us assume that we believe that the state of our multi- or bipartite quantum 
system is an entangled state and let us assume that we are interested in 
detecting the entanglement of over some bipartite A-B cut Ha <8> TIb- For 
pure states, there always exists a witness H of the form 

H = aP + (l-a)(l®T)(Q), (35) 

where P > and Q > 0, see Section 2.2. An optimal witness (see Theorem 
2, Ref. [LKCH]) in this class has a = and has the property that the 
operator Q has no product states in its range. To optimize with respect to 
the state we choose Q = iVVwn) (VVmiJ where |VVwn) is the eigenvector of 
(1 (g) T)(\ip)(ij}\) which has the smallest eigenvalue (which is negative). Such a 
witness is optimal in the sense of Ref. [LKCH] . The optimality of this choice 
with respect to follows from 

Tr(l®T)(Q)^)(V| = TrQ(l®T)(|V><V|) > A min , (36) 
since Q > and TyQ = 1 due to normalization. 

To see the explicit form of such a witness, we write in the Schmidt decom- 
position = J2i y/Kl^i) ® \h) and take the partial transpose in a fixed basis 
|0),|1),.... We have 

(l^TX^y^KKa,-! ® l&^-l) =J2^X j \a i )(a j \ <g> (37) 

which has eigenvectors and corresponding eigenvalues {|cij, &*), Aj} and for z ^ 

3, 

^{\a i )®\b^±\a j )®\b*)), iV^}- (38) 

We choose the eigenvector with the most negative of the eigenvalues, let us 
call it Umin — ~ maxj^j <J XiXj. The optimal ?/>- witness is equal to 

H* = ^{\ai,bj)(ai,bj\ + \aj,bi)(aj,bi\ - \ai,bi){aj,bj\ - |o,-, bj)(ai, &i|).(39) 



2 We will only be optimizing amongst the decomposable witnesses to keep things 
as simple as possible. 



19 



for the pair corresponding to p m in- Let us take a simple example. 
Example 5 For the state = cos#|00) + sin#|ll), the optimal witness is 

H* = ±(|0, 1)(0, 1| + |1,0)(1,0| - |0,0><1, 1| - |1, 1)(0,0|), (40) 
and fimin = — Vcos 9 sin 6 . 

For bipartite quantum systems consisting of 2 qubits or 1 qutrit + 1 qubit, the 
entanglement witness is always decomposable, i.e. of the form H = aP + (1 — 
a)(l <g>T)(Q) (see Section 2.2). As for pure states, the witness is optimal when 
a = and Q has no product states in its range. To optimize for p among such 
witnesses, we find the eigenvector IVVmm) [f] of (1 <g> T)(p) with the smallest 
eigenvalue p min and we choose 

#* = (l®T)(|Vw)<VwJ)- (41) 



4-2 Mixed State Entanglement in Higher Dimensions 

Let p be a (multipartite) mixed state in dimensions more than H, 2 (£> H.2 or 
® 7~^3 whose entanglement we wish to detect over a bipartite cut A-B. 
If p violates the Peres-Horodecki criterion then the methods in the previous 
Section can be applied to determine an optimal entanglement witness. When 
p has the PPT property and is believed to be entangled, then p has bound 
entanglement and we will need an indecomposable entanglement witness for p 
(see Section 2.2). In Ref. [LKCH] a method was developed to optimize a given 
indecomposable entanglement witness. The problem is a lot harder than for 
decomposable witnesses: a generic form for an optimal witness is not known. 
Nonetheless we can sketch a procedure for finding an optimal indecomposable 
witness which detects p: 

(1) Find a decomposable witness H, Eq. (10), for which Tr Hp = with 
TrH = 1. 

(2) Choose as a starting point the indecomposable witness H' = (if— el) / (1 — 
ed) where e = inf^^^i, ^ 2 | H \ipi, ^2) which ensures that H' is a (nor- 
malized) witness. Here d is the total dimension of the quantum system. 

(3) Optimize H' to H* with the methods in Ref. [LKCH], taking into account 
the optimality with respect to p, see Definition 4. The optimized witness 



3 This eigenvector is always entangled, since for all product states I^a^b), 
fe| (1 T)(p) fe) = Tr p (1 ® T)(|^, > 0. 
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relates to if' as 



if, = (H'-\ m D m )/(l-\ m ), 



(42) 



where A* is a constant depending on .D* and if' (see Eq. (13) in Ref. 
[LKCH]) and where -D* is a normalized decomposable witness, Eq. (35), 
with the property that for all product states \ipi,ip2) such that (1) 
(■01, V2I H' |^i, ^2) = we have (01, V^l -D* l^i, ^2) = 0, (2) A* and 
are such that if* is an optimal indecomposable entanglement witness (no 
more decomposable witnesses can be subtracted from it) and (3) when 
under constraints (1) and (2) there is still freedom of choice in A* and 
f}*, we choose an ff* which is optimal with respect to p, i.e. Tr H*p is 
minimal. 

It is not guaranteed that this procedure will lead to a p-optimal indecompos- 
able witness according to Definition 4. The method does however ensure that 
the witness is optimal as well as detecting the entanglement of p. 

4-3 Multipartite (Bound) Entanglement 

The entanglement witness framework carries over to multipartite states. This 
framework is particularly useful when entanglement is to be detected in a 
bound entangled state such as the first example pshift S) Eq- (32), in Section 
3.4 which is separable over all bipartitions. For this state an entanglement 
witness was found in Ref. [HHH], in analogy to the construction in Ref. [Ter]. 
Even though various Bell inequalities and entanglement witnesses exist for 
multipartite entanglement (Sections 3.1 and 3.2), a more elaborate theory as 
in the bipartite case is still lacking. We would also like to refer the reader to 
Ref. [DC] for alternative ideas on the experimental detection of entanglement 
in multiqubit states. 



4-4 The vicinity of the trial state 

Often the quantum system under consideration will not exactly be in the 
desired state p. It is therefore essential that the entanglement witness that we 
chose is robust, in the sense that it detects as much entanglement as possible 
in the neighborhood of p (provided that the neighborhood is entangled). This 
property is guaranteed in the following manner. Let p' be all states in the 
'vicinity' of p, namely for a given e, A = p' — p 



Alii < e, 



(43) 
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where ||.||i is the tracenorm, i.e. ||A||i = TrV A^A. Letting if* be the optimal 
p-witness, we have 

Tr iV = Tr H+p + Tr if, A. (44) 

The last quantity can be bounded, using the Schwarz inequality |Tr A*B\ < 
VTyATAVTiBW and also VTr A^A < Tr VA^A. This gives 

Tr H+p - e ^TrH* < Tr H*p' < Tr H*p + e ^TrH*. (45) 

For the optimal decomposable entanglement witnesses given in Section 4.1 we 
have ^JTtH^ = 1, irrespective of the dimension of the system. Thus when 
p' is close to the state p, the optimal witness if* for p will also detect the 
entanglement in p'. We will now consider a specific example in which all states 
in a certain class are detected by one entanglement witness: 

Example 6 Let be any maximally entangled pure state in a bipartite space 
of total dimension d. Instead of the physical state of our quantum system 
is 

Pp =p\m)(m\ + l -^Pi. (46) 

As long as (l<g>T)(p p ) ^ ; the optimal^ -witness H* = {^-®T)(\ip ^ min ) {ip ilmin \) 
will detect the entanglement of p p . The eigenvector with the minimal eigenvalue 
of (1 <E> T)(p p ) is the same for all p G [0, 1]. Furthermore the density matrix 
p p has the property that it is separable if and only if it satisfies the Peres- 
Horodecki criterion [HH99]. Therefore the entanglement in all (entangled) p p 
is witnessed by H*, since Tr H^pp is negative as long as p p violates the Peres- 
Horodecki criterion. 



4-5 Measurement of the witness 

In principle, the entanglement witness method has the advantage that only one 
observable, the entanglement witness, needs to be measured. In practice, the 
measurement of this observable may be done by a series of local measurements. 
Consider for example the entanglement witness of Example 5, in terms of the 
Pauli-matrices it reads 

if. = ^{I®I + Y ®Y -X®X- Z®Z). (47) 
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A measurement of all the Pauli matrices on both sides is needed to measure this 
observable by local measurements. At this point the advantage over basic state 
tomography becomes somewhat questionable. The entanglement witness will 
in general be a nonlocal observable. But if we allow for nonlocal measurements, 
or rotate the state by a nonlocal rotation U prior to testing such that URJJ^ 
is local, there exists a real advantage over state tomography. A CNOT gate 
where the first qubit is the control and the second is the target, will rotate if* 
in Eq. (47) to a product of operators 

CNOT^ B H* CNOTJ^b = -(1 - X) <g> (1 - Z)). (48) 

Thus some quantum computation power is needed to measure the entangle- 
ment witness efficiently For large systems one may ask whether such quantum 
computation can be performed efficiently, in polynomial time in the number 
of qubits, for states which can be built efficiently. With this open question we 
will conclude our review. 
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A Indecomposable Positive Linear Maps and Hilbert's 17th prob- 
lem 

In Ref. [Cho75] it was shown how certain indecomposable positive linear maps 
present answers (in the negative) to Hilbert's 17th problem. Hilbert's 17th 
problem asks whether all positive semidefinite homogeneous polynomials are 
sums of squares of homogeneous polynomials (see Ref. [Rez] for some history 
on the problem). Indecomposable positive maps with real coefficients present 
counterexamples for polynomials which are real biquadratic forms. The con- 
struction is the following. 

Given is a positive indecomposable map C: B(Ti n ) — > B(7i m ) which is known 
not to be completely positive. Furthermore, the coefficients of C in some fixed 
basis are real and therefore C maps real (symmetric) matrices in B(7i n ) 
onto real symmetric matrices in B{TL m ). Note that matrix transposition acts 
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as the identity on real symmetric matrices and therefore we restrict ourselves 
to indecomposable maps. 

We define a positive semidefinite symmetric biquadratic form in the following 

way: 

F(x 1 ,...,x n ;y 1 ,...y m ) = (y\ (C(\x)(x\) \y), (A.l) 

where \x) and \y) are unnormalized vectors with coefficients \x) = J2i x i\i) an d 
\y) — Y^iVi\i)j Xi EH and G R. The positivity of the map £ guarantees the 
positive semidefiniteness of F(x±, . . . , x n ; yi, . . . y m ) for all x±, . . . , x n , y±, . . . , y m 
or 

F(x u ...,x n ;y 1 ,...y m )= J2 £ij,ki x i x jVkyi > 0. (A.2) 

i y j,k,l 

Now we pose Hilbert's question: is every positive semidefinite symmetric bi- 
quadratic form a sum of squares, i.e. 

F(x u ...,x n ;y 1 ,...y m ) = ^2(G t (xi, . . . , x n ; y u . . . y m )) 2 , (A.3) 

t 

where G t (xi, . . . , x n ; yi, ■ ■ ■ y m ) is a symmetric bilinear form, i.e. 

G t (x 1 ,...,x n ;y 1 ,...y m ) = Y^g\jXiyj. (A.4) 

Assume that the equality in Eq. (A.3) holds. We define a set of real operation 
elements A t with 

{j\AM=gh, (A.5) 

so that 

J2(Gt(xi, . . . ,x n ;y±, . . . y m )f = J2iv\ M x ) (A^W)- ( A - 6 ) 
t t 

This would imply that the map C has a decomposition in terms of the oper- 
ation elements A t and therefore £ is a completely positive map. Conversely, 
when a map C is completely positive, the corresponding positive semidefinite 
symmetric biquadratic form is a sum of squares. Instead of the positive inde- 
composable map C, we could have defined the biquadratic form in terms of 
an (indecomposable) entanglement witness H with real coefficients: 

F(xi, . . . , x n ; y u . . . , y m ) = (y, x\ H \y, x) > 0, (A.7) 
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for all real vectors \y,x). It was shown how to construct (real- valued) entan- 
glement witnesses for every (real) unextendible product bases (UPB) in Ref. 
[Ter]; furthermore the graph-theoretic UPB construction in Ref. [AL] always 
has a realization with real vectors. 
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